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Abstract 

In contrast to the canonically conjugate variates q,p representing 
the position and momentum of a particle in the phase space distri- 
butions, the three Cartesian components, J x ,J y , J z of a spin-j sys- 
tem constitute the mutually non-commuting variates in the quasi- 
probabilistic spin distributions. It can be shown that a univariate 
spin distribution is never squeezed and one needs to look into either 
bivariate or trivariate distributions for signatures of squeezing. Sev- 
eral such distributions result if one considers different characteristic 
functions or moments based on various correspondence rules. As an 
example, discrete probability distribution for an arbitrary spin-1 as- 
sembly is constructed using Wigner-Weyl and Margenau-Hill corre- 
spondence rules. It is also shown that a trivariate spin-1 assembly 
resulting from the exposure of nucleus with non-zero quadrupole mo- 
ment to combined electric quadrupole field and dipole magnetic field 
exhibits squeezing in cerain cases. 
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1 Introduction 

A spin-j' system can be represented by an n x n density matrix where n = 
2j + 1. Any spin-j density matrix can be parametrized jTj in terms of n 2 — 1 
independent, spherical tensor parameters t k as 

2j +k 



Z J 1 k=0q=-k 

where 



t k = Tr{pr k q {J)). (2) 
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Spherical tensor operators tHj) can be obtained by Weyl's construction |2] 
using spin operators J x , J y , J z as 

r k q {J) = N k ^J.V) k r k Y^e^) (3) 

where are the normalization factors given by 



ok 



47r(2j-fc)!(2j + l) 
k\\ (2j + k + l)\ { ] 

which are in agreement with the earlier usage of Ramachandran and Umergee [3] 
and with the Madison convention [I]. Here Y q (6, 0) are the well known spher- 
ical harmonics. The spherical tensor operators satisfy the orthogonality re- 
lations, 

Tr(T*(J)T*\?))=n5 kk , qq/ (5) 

where 

rf(J) = (-l)V g (J). (6) 
Hermiticity of the density matrix together with Eq.(6) demands that 

The set of n 2 — 1 spherical tensor operators T q (J) together with the identity 
matrix, form an orthonormal basis of the linear vector space of complex nxn 
matrices that act on the n dimensional spin space. Thus p may be completely 
determined by finding the expectation values of these n 2 — 1 operators in the 
state p. 

Consider the example of a spin-1 density matrix. It can be completely 
characterized by the spherical tensor parameters ij, t^, t^, t\ x , tj_ 2 and the 
corresponding spherical tensor operators, r q and r 2 can be written as 

r 1 q (J) = N 1 - j \j.V)rY q \9A) (8) 
r 2 q (J) = N 2 /(J.V) 2 r 2 Y 2 (e, < j ) ). (9) 

" ■ J 4 ■ J i n \l' ] (10) 



and 

Explicitly, 
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<=^ + 4 + W- 1 ) (12) 

= TN 2 ~/(J x ^ + J y ^ + J^f^[z{x ± iy)] (13) 
'9y 



and 



*5= N ^\fI Tr P J z (15) 

4 1= TN^y/^Tr(p{J x ±iJ Y )) (16) 

*o = iV^V^T^SJl-jO' + l))) (17) 

4i= TNy 1 y/£Tr(p((J z J x + J x J z )±i(J z J Y + J Y J Z ))) (18) 
*±2 = N^^Tr(p((J x - J Y ) ± 2*(J X J y + J y J x ))). (19) 

Here one could easily identify as linear combinations of I s * and 2 nd 

order Cartesian moments of the operators J x ,Jy,Jz- Further, computing all 
the spherical tensor parameters characterising a spin-j density matrix, one 
could see that all the t£'s can be constructed by taking appropriate linear 
combinations of the k th order cartesian mixed moments p abc , (a + b + c = k) 
of the operators J x , Jy, Jz provided the mixed moments are defined by the 
relation 

^ = {a ^ C l cy Tr[pS ww (J a xJYJ c z)} (20) 

where S ww () denotes the the complete symmetrizer, also known as Wigner- 
Weyl symmetrizer, with terms. In classical probability theory mixed 

moments are defined through characteristic function $(/) as 

aa ab nc 

» abc = (-*r b+c WaWWc Hi)L=i y ^, (2i) 

w x y z 
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Therefore one can identify the characteristic function leading to the multi- 
variate moments in Eq.(20) to be 

<$> ww (I)=Tr(pe lTJ ). (22) 

Note that this function can be obtained by the Wigner-Weyl prescription, 

e ii-x e ii-j ( 23 ) 

or eqivalently, 

X a Y b Z c = (Q + fc + c)! S ^( J x4^)- (24) 

Here X, Y, Z are the classical variates corresponding to the quantum me- 
chanical operators Jx,Jy,Jz respectively. Taking the fourier transform of 
the characteristic function, <& WW (I) one can obtain the Wigner-Weyl proba- 
bility density function as 

fUX)= Poo Poo Poo^ I $_(/)e-^ (25) 

= j^EtoEU-^K^ k q (X)l k (X) (26) 

where 



* W = ^ t c ^r, £ di rV*-> (27) 



and V k q (X) are the spherical components of the gradient operator and C(jkj; mOm) 
are the well known Clebsch-Gordon coeffecients. One could easily verify that 

J °^ = T^^ s{x - m) (28) 

and 

= ^9S? c( *' imom) ^ i(A ' - mri<x - m)l (29) 

for k > 1. Here S(X — m) denotes the Dirac 5— function and e(X — m) 
denotes the step function 

1 poo piat 

e(a) = — di—. (30) 
2m J-oo t 



Here e(a) = 1 if a > and e(a) = if a < 0. 

Thus the problem of finding a quantum mechanical density function re- 
duces to finding the quantum mechanical equivalent of the classical char- 
acteristic function or the classical moments. Owing to the noncommutative 
nature of the spin components Jx, Jy and Jz, one can realize several quantum 
mechanical definitions of the characteristic function $(/) or mixed moments 
fi abc based on different operator orderings. Once the moments are obtained, 
instead of a continuous probability distribution function, one can obtain a 
discrete probability mass function. The procedure for obtaining a discrete 
probability distribution for a spin-j system is outlined in Sec. 2 . Rules asso- 
ciating classical quantities to quantum mechancial operators are called corre- 
spondence rules. There are many such rules|6j. Of these, the Margenau-Hill 
and the Wigner-Weyl procedures are well studied in the discussion of quan- 
tum phase distributions and quantum spin distributions jZj,jS]. As an exam- 
ple we would like to compute the probability mass function for any general 
spin-1 system using both Margenau-Hill and Wigner-Weyl correspondence 
rules in Sec. 3. In Sec. 4 we discuss the squeezing of univariate and multivari- 
ate spin distributions. 

2 Discrete probability distributions of spin-j 

system 

Margenau-Hill correspondence rule is given by 



X Y Z — > , [Jx^y^z + Jx^z^y + Jy^x^z + Jy^z^x (32) 




(31) 



or 



I jc ja jb I jc jb ja 1 
^ J Z J X J Y + J Z J Y J X\- 



Therefore the mixed moments in this scheme are given by 



abc 



1 

3! 



Tr[pS MH {J x J b Y J c z }} 



(33) 
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where Smh{} stands for the Margenau-Hill symmetrizer as given on the right 
hand side of Eq.(32). Therefore, the Marganeau-Hill characteristic function 
can be written as, 

<& MH = ^Tr(p[e iIxJx e iI y jY e iIzJz + e iIxJx e ihJz e iI y jY + (34) 

gilyjy gilxJx gilzJz _|_ gilyJy gilzJz gilxJx _|_ qHzJz qUxJx gUyJy 
j^gilzJz gUyJy gilxJx^ 

Similarly, the Wigner-Weyl correspondence rule is given by Eq.(23) or 
equivalently by Eq.(24) and the characteristic function by Eq.(22). Thus the 
mixed moments in this scheme can be computed using Eq.(20). 

Note that the classical variates XY Z have the same diagonal form as 

J x ^ '~^y i ^ z ; theft 



V 



+3 



U - i) 



-u - i) 



\ 



and the space of these random vairables 



has (2j + l) 3 points. Because of the algebra of spin matrices there are only 
(2j + l) 2 — 1 independent moments and all the others can be expressed in 
terms of these moments only. 

If P{m x m y m z ) is the probability mass function, then by definition the 
multivariate moments are given by, 



^(xyz) = E E E 



m x m y m z ) 



(35) 



ra, m y m z 



Here, a, f3, 7 = 0, 1, 2, ....2j and m x , m y , m z are the values taken by the ran- 
dom variables X, Y, Z respectively. Since X, Y, Z are the classical random 
variables corresponding to the quantum mechanical operators Jx, Jy , Jz, for 
a spin-j system, m x , m y , m z = +j, , —j. Eq.(35) can be written as 



and explicitly , 



^(XYZ)=J2«, 
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m,x 



(36) 



,107 



,,2307 



1 1 

i 3-1 
J 2 U - 1) 2 



1 1 

-a - 1) -3 
(-i) 2 



/i = VA 



{-if 



3 — 1 



(37) 
(38) 
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V is the well known Vandermonde matrix. This equation can be inverted to 
get A^ x in terms of the known moments // 0/37 , /i 1/37 , •, •, /i 2jf/37 . But 

"^mL = m y^m x rn y (^9) 
m y 

where 

S^m, = 5I(m z ) 7 P(m :r m 3/ m z ). (40) 

Thus P{m x m y m z ) can be computed easily by successive inversions. Observe 
that different sets of moments yield different probability distributions. 



3 Trivariate probability mass function for spin- 

1 system 

In the case of spin-1 matrices, it has been shown by Weaver [S] that a par- 
ticularly useful set of algebraic relations arises when one defines the set of 
matrices, S l5 S 2 , S 3 to be 

Si = J!~ J'y (41) 

XI2 Jx,Jy H~ JyJx 

s 3 = J z . 

Using the relation 

JiJjJk Jk'JjJi — $ijJk $jkJi (42) 

the algebra of Ej's has been shown to be 



y^2 v^2 v^2 
Zjj — Z^ 2 — Z^ 3 — 




(43) 















which combine to give 



EjSj = + ie iifc S fc (44) 



exactly same as the Pauli spin matrix algebra. Here represents the Levi- 
Civita symbol. 
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Since X,Y,Z are given by the diagonal matrix ^ « i ^ , the space of 

our random variables have only 27 points and we need exactly 27 moments 
// jfe ; k = 0, 1, 2 to compute the probability mass function. 
Using the relations (41) to (44) extensively together with 

Tr[p{Jl + Jy + Jl)} = 2, (45) 

it can be shown that, of the 27 moments, only 8 are independent. We choose 
the independent moments to be p 001 , p ow , p wo , p 200 p 020 , p wl , p no , p 011 . The 
remaining 19 moments can be computed interms of the above 8 independent 
moments. 

Thus the 27 points of the probability mass function both in Margenau-Hill 
and Wigner-Weyl prescriptions can be obtained interms of the 8 independent 
moments as 

1 1 

PMH(m x m y m z ) = — ^ m i Tr[ / 9j i ] + — ^ m i m k Tr[p{3 i 3 iL + J k Ji)], 

i=x,y,z k=ii=x,y,z 

(46) 

when all the 3 arguments are non-zero 

P M H{m x m y m z ) = ± ^i=x, y ,z miTr[p3^\ + \[E i=x ,y,z m%Tr\pJf] - 1] (47) 

+ ik mim k Tr[p{J iJ k + J k J ; )], 

when two of the arguments are non-zero, 

PMH{m x m y m z ) = — ^ m;Tr[pJi] (48) 

i=x,y,z 

when two of the arguments are zero and 

PMH(m x m y m z ) = (49) 

where m x = m y = m z = 0. 

Wigner-Weyl probability mass function turns out to be 

Pww(m x m y m z ) = -!- mTr[pJi\ + — ^ m;m fc Tr[p(JiJ k + J k Ji)] 



120 . ^ Lr J 96 , 

i=x,y,z k^=i=x,y,z 



(50) 
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when all the 3 arguments are non-zero, 

Pww(m x m y m z ) = ^- «^r[pJi] + — ^ m;m fc Tr[p(JiJ k + J k Ji)] 

i=x,y,z kj^i=x,y,z 

(51) 

when two of the 3 arguments are non-zero, 

15 1 

Pww(m x m y m z ) = - J2 m i Tr [pJi] + J2 m i Tr \p(3i 2 )] - g ( 52 ) 

i=x,y,z kj^i=x,y,z 

when two of the 3 arguments are zero and 

P W w{m x m y m z ) = — (53) 

where m x = m y = m z = 0. 
Notice that 

PMH{m x m y m z ) = 1. (54) 

m x ,m y ,m z 

and 

P W w{m x m y m z ) = 1. (55) 

mx,m v ,mz 

From the expressions for the probability mass function it is evident that 
they can assume negative values also. Summing over any one of the indices 
in the expressions for probability distributions one can obtain the biavriate 
probability mass function which is capable of describing the aligned, biax- 
ial spin-1 system which can be generated by exposing a spin-1 nucei with 
non-zero electric quadrupole moment to external electric quadrupole field 
generated by a suitable crystal lattice [Hj. Summing over any of the two 
indices we can get the univariate probability distribution. As expected both 
Margenau-Hill and Wigner-Weyl distributions lead to the same univariate 
distribution. Such distributions of any spin-j system represent oriented spin 
systems jTUj. For instance an oriented assembly reults when a spin-j system 
is exposed to external dipole magnetic field. Such probability distributions 
are semi positive definite and are nearest to the classical probability distri- 
butions. Where as , the bivariate and trivariate distributions are generally 
quasi probabilistic and the corresponding spin assemblies are called non- 
oriented|IT . Thus any spin-j system has to belong to one of the mutually 
exclusive classes viz., oriented or non-oriented systems. It is interesting to 
study the squeezing property of the univariate, bivariate and trivariate spin-j 
systems. 
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4 Spin squeezing 



One may consider the spin system to be squeezed in J x or J y provided 

AJ 2 x (\\TrpJ z \ (56) 

or 

AJ 2 y (±\TrpJ z \ (57) 

where A Jy and A Jy represent the variances in J x and J y respectively It has 
already been noted by Kitagawa and Ueda[12j that this definition implies 
that a coherent spin state is squeezed if placed in an appropriate coordinate 
system. To overcome this difficulty, they consider a spin-j system to be made 
up of 2j elementary spinors and choose a coordinate system whose z-direction 
coincides with the mean spin direction. 

For a mixed system, the vector polarization is defined as 

Kitagawa and Ueda l2. refer to this as the mean spin direction. 

The spherical components of vector polarization P are related to through 



P. = ^^±% (59) 

Choosing the Z-axis along P one is left with infinite choices of X and Y axes. 
Now make a rotation about P such that the covariance 

{AJ x J y f = Tr (p( J ^ + W) _ Tr (pJ x )Tr(pJ y ) = 0. (60) 
In this resulting frame, the Schrodinger inequality 

(Aj x nAj y n-Aj x j y r > gM^w (61) 

and the Heisenberg- Robertson inequality [TH]. 

{AJ x f{AJ y f > (62) 
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are of the same form and the definition of spin squeezing is unambiguous ie., 



{AJ X 



< 



or 



(AJ y ) 2 < 



\Tr(pJ z 
2 

\Tr(pJ z . 



In terms of spherical tensor notation, the squeezing criterion becomes 



1 + 



3(2j + 3)(2j- 
j(j + 1)40 



11/2 




-ti) < - 



11/2 



j(j + 1) 



(63) 
(64) 

(65) 



This condition is never satisfied by oriented systems [T4*]. in other words uni- 
variate systems. One has to consider bivariate or trivariate systems. The 
simplest bivariate system |T3] one could think of is a spin-j assembly of nuclei 
with non-zero electric quadrupole moment exposed to external electric and 
quadrupole field generated by electrons in suitable crystal lattice site. Here, 
£q = and the squeezing criterion is never satisfied. Therefore one needs to 
consider a trivariate system. If the above bivariate system is exposed to an 
external magnetic field along the z-axis of the Principle Axes Frame(PAF) , 
the resultant spin assembly will need all the three variates J x , J y , J z for its 
description^]- In the case of spin-1 assembly , there are situations when 
the assembly exhibits squeezing and a few values are presented in table. 1. 

Table 1: Trivariate, squeezed spin-1 system specified by its 
non-zero spherical tensor parameters 



tl 


t 2 


t 2 


(AJ X ) 2 


\Tr(pJ z )\ 


-0.7506 


0.495 


-0.4453 


0.2929 


0.3064 


-0.6298 


0.4737 


-0.5307 


0.2486 


0.2571 


-0.7506 


0.4526 


-0.4453 


0.3028 


0.3064 



5 CONCLUSIONS 

In this paper it has been shown that any general spin-j assembly can be char- 
acterised by a trivariate probability distribution function with J x ,J y ,J z as 
the variates. The distribution needs n 2 — 1 independent moments of J x , J y , J z 
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for its complete description. Because of non-commutativity of these variates 
there exist various correspondence rules which associate different quantum 
mechanical expressions to the same classical expression for moments. Using 
two such correspondence rules viz, Wigner-Weyl and Margenau-Hill rules, 
we have computed the probability mass function to arbitrary spin-1 systems. 
These distributions are in general found to be quasi-probabilistic. Starting 
from moments we have also outlined the method to obtain the probability 
mass function by the method of inversion. Since univariate distributions cor- 
respond to oriented spin assemblies, they are never squeezed. Thus one has 
to look for signatures squeezing in the bivariate or trivariate distributions. 
In the case of spin-1 assembly, the trivariate distribution which describes a 
spin-1 nucleus like H 2 , N 14 exposed to combined electric and magnetic field 
shows squeezing in certain cases. 
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